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^ ! Abstract 

\ By introducing a concept of thermal expansion (TE) of a Josephson junction as an elastic 

' response to an effective stress field, we study (both analytically and numerically) the temperature 

o ' 

+-> 



and magnetic field dependence of TE coefficient a in a single small junction and in a square array. 
In particular, we found that in addition to field oscillations due to Fraunhofer-like dependence 



of the critical current, a of a small single junction also exhibits strong flux driven temperature 
oscillations near Tc- We also numerically simulated stress induced response of a closed loop with 
finite self-inductance (a prototype of an array) and found that a of a 5 x 5 array may still exhibit 
temperature oscillations provided the applied magnetic field is strong enough to compensate for 



the screening induced effects. 

PACS numbers: 74.50. +r, 74.62.Fj, 81.40.Jj 
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I. INTRODUCTION 



Inspired by new possibilities offered by the cutting-edge nanotechnologies, the experimen- 
tal and theoretical physics of increasingly sophisticated mesoscopic quantum devices (heavily 



Decoming one of the most exciting and 
JfiQ farther .views charge 



based on Josephson junctions and their arrays) is 
rapidly growing areas of modern science (see, e.g., [l 
and spin effects in mesoscopic 2D Josephson junctions and quantum-state engineering with 
Josephson-j unction devices). In particular, a remarkable increase of the measurements tech- 
nique resolution made it possible to experimentally detect such interesting phenomena as 
flux avalanches {s| and geometric quantization [6] as well as flux dominated behavior of heat 
capacity [7] in Josephson junctions (JJs) and their arrays (JJAs). 

At the same time, given a rather specific magnetostrictive [sj and piezomagnetic re- 
sponse of Josephson systems, one can expect some nontrivial behavior of the thermal expan- 
sion (TE) coefficient in JJs as well. Of special interest are the properties of TE in applied 
magnetic field. For example, some superconductors like Bai- x K x BiOs, BaPb x Bi\_ x O?, 
and La 2 - x Sr x CuO \ were found [h]] to exhibit anomalous temperature behavior of both 
magnetostriction and TE which were attributed to the field-induced suppression of the su- 
perstructural ordering in the oxygen sublattices of these systems. 

By introducing a concept of TE of Josephson contact (as an elastic response of JJ to 
an effective stress field), in the present paper we consider the temperature and magnetic 
field dependence of TE coefficient a(T, H) in a small single J J and in a single plaquette (a 
prototype of the simplest JJA). In a short contact, the field-induced a(T,H) is found to 
exhibit strong temperature oscillations near Tq- At the same time, in an array (described 
via a closed loop with finite self-inductance) for these oscillations to manifest themselves, 
the applied field should be strong enough to overcome the screening induced self-field effects. 

II. THERMAL EXPANSION OF A SMALL JOSEPHSON CONTACT 

Since thermal expansion coefficient a(T,H) is usually measured using mechanical 



dilatometers 



111 ] . it is natural to introduce TE as an elastic response of the Josephson 

Namely, we define the TE coefficient (TEC) 



contact to an effective stress field a 
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at(T, H) as follows: 

«T,B)=*L (1) 

where an appropriate strain field e in the contact area is related to the Josephson energy Ej 
as follows (V is the volume of the sample): 

~dEj 



1 
V 



da 



(2) 

<7=0 



For simplicity and to avoid self-field effects, we start with a small Josephson contact of length 
w < \j (Aj = ^J&o/{iodj c is the Josephson penetration depth) placed in a strong enough 
magnetic field (which is applied normally to the contact area) such that H > $o/27rA jd, 
where d = 2A^ + 1, Xl is the London penetration depth, and t is an insulator thickness. 

The Josephson energy of such a contact in applied magnetic field is governed by a 
Fraunhofer-like dependence of the critical current [13J: 

/ simp \ 
Ej = J I 1 — cos^ol, (3) 

where if = 7r$/$o is the frustration parameter with $ = Hwd being the flux through the 
contact area, cpo is the initial phase difference through the contact, and J oc e~*'^ is the 
zero-field tunneling Josephson energy with £ being a characteristic (decaying) length and t 
the thickness of the insulating layer. The neglected here self-field effects (screening) will be 
treated in the next Section for an array. 

Notice that in non-zero applied magnetic field H, there are two stress-induced contri- 
butions to the Josephson energy Ej, both related to decreasing of the insulator thickness 
under pressure. Indeed, according to the experimental data [l2j, the tunneling dominated 
critical current I c in granular high- Tc superconductors was found to exponentially increase 
under compressive stress, viz. J c (cr) = I c (0)e Kcr . More specifically, the critical current at 
a = 9kbar was found to be three times higher its value at a = 1.5kbar, clearly indicating 
a weak-links-mediated origin of the phenomenon. Hence, for small enough a we can safely 
assume that [9] t(a) ~ t(0) (1 — (3a/ ao) with ao being some characteristic value (the param- 
eter (3 is related to the so-called ultimate stress a m as (3 = a /a m ). As a result, we have the 
following two stress-induced effects in Josephson contacts: 

(I) amplitude modulation leading to the explicit stress dependence of the zero-field energy 

J(T, a) = J(T, 0)e^ /ao (4) 



with 7 = pt(0)/£, and 



(II) phase modulation leading to the explicit stress dependence of the flux 



$(T,H,a) = Hwd(T,a) 



(5) 



with 



d(T, a) = 2\ L {T) + t(0)(l - (3a/a ) (6) 

Finally, in view of Eqs.(l)-(6), the temperature and field dependence of the small single 
junction TEC reads (the initial phase difference is conveniently fixed at <po = n): 

dF(T, H) 



a(T, H) = a(T, 0) [1 + F(T, H)} + e(T, 0)- 



dT 



where 



with 



F(T,H) 



Sill Lf 



<P{T,H) 



<f d(T, 0) \ if 
ir$(T,H,0) H 



COS If 



a(T, 0) 



$ H (TY 

de(T, 0) 
~ dT ' 



and 



(7) 

(8) 

(9) 
(10) 

(11) 



Here, H (T) = $ /vrwc/(T, 0) with d(T, 0) = 2A L (T) + t(0). For the explicit temperature 
dependence of J(T, 0) = QqIc(T) /2tt we use the well-known [14| analytical approximation 



of the BCS gap parameter (valid for all temperatures), A(T) = A(0) tanh \ 2.2\J Tc ~ T 
with A(0) = 1.76/cb?c which governs the temperature dependence of the Josephson critical 
current 



I C {T) = I c (0) 



A(T) 
A(0) 



tanh 



A(T) 
2ksT 



(12) 



while the temperature dependence of the London penetration depth is governed by the 
two-fluid model: 

\l{T) = . Al(0) (13) 

From the very structure of Eqs.(l)-(9) it is obvious that TEC of a single contact will 
exhibit field oscillations imposed by the Fraunhofer dependence of the critical current Ic- 
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FIG. 1: (Color online) Temperature dependence of the normalized flux driven strain field 
e(T, /)/e(0, 0) in a single short contact for different values of the frustration parameter / = H/Hq(Q) 
according to Eqs.(l)-(13). 

Much less obvious is its temperature dependence. Indeed, Fig. [1] presents the temperature 
behavior of the contact area strain field e(T, /) (with i(0)/£ = 1, £/Al(0) = 0.02 and (3 = 0.1) 
for different values of the frustration parameter / = H/H Q (0). Notice characteristic flux 
driven temperature oscillations near Tq which are better seen on a semi-log plot shown in 
Fig. [2] which depicts the dependence of the properly normalized field-induced TEC a(T, f) 
as a function of 1 — T/Tq for the same set of parameters. 

III. THERMAL EXPANSION IN THE PRESENCE OF SCREENING CURRENTS 

To answer an important question how the neglected in the previous Section screening 
effects will affect the above-predicted oscillating behavior of the field-induced TEC, let us 
consider a more realistic situation with a junction embedded into an array (rather than an 
isolated contact) which is realized in artificially prepared arrays using photolithographic tech- 
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FIG. 2: (Color online) Temperature dependence of flux driven normalized TEC a(T, /)Tc/|e(0, 0)| 
in a single small contact for different values of the frustration parameter / = H/Hq(0) (for the 
same set of parameters as in Fig.l) according to Eqs.(l)-(13). 



nique (that nowadays allow for controlled manipulations of the junctions parameters 15]) 



Besides, this is also a good approximation for a granular superconductor (if we consider it 



as a network of superconducting islands connected with each other via Josephson links [13]). 
Our goal is to model and simulate the elastic response of such systems to an effective stress 
a (described in the previous Section for an isolated contact). For simplicity, we will con- 
sider an array with a regular topology and uniform parameters (such approximation already 
proved useful for describing high-quality artificially prepared structures |6|]). 

A. Model equations for a planar square array 

Let us consider a planar square array as shown in Fig. [3J The total current includes 
the bias current (flowing through the vertical junctions) and the induced screening currents 
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FIG. 3: Sketch of an array. The junctions perpendicular (parallel) to the bias are called horizontal 
(vertical), (a) The node is shown as a circle in the left bottom corner of a plaquette; (b) a 
single plaquette (the elementary unit of the circuit) along with the circulating current; and (c) the 
lumped elements circuit for a small junction. 

(circulating in the plaquette [3]). This situation corresponds to the inclusion of screening 
currents only into the nearest neighbors, neglecting thus the mutual inductance terms 
Therefore, the equation for the vertical contacts will read (horizontal and vertical junctions 
are denoted by superscripts h and v, respectively): 

^ + 2^r +I ^ = AI '^ h ' (14) 

where Alf- = If- — If , • and the screening currents I s obey the fluxoid conservation con- 
dition: 

- 4>h + - ti d + Aw = 2 -^7 - (15) 
Recall that the total flux has two components (an external contribution and the con- 
tribution due to the screening currents in the closed loop) and it is equal to the sum of 
the phase differences describing the array. It is important to underline that the external 
flux in Eq. (|15|) . rj = 2nQ ext /$ 0; is related to the frustration of the whole array (i.e., this is 



the flux across the void of the network 



18 



191]). and it should be distinguished from the 



previously introduced applied magnetic field H across the junction barrier which is related 
to the frustration of a single contact / = 2-nHdw /<3> an d which only modulates the critical 



current Ic(T, H, a) of a single junction while inducing a negligible flux into the void area of 
the array. 



B. Stress induced effects 

For simplicity, in what follows we will consider the elastic effects due to a uniform (ho- 
mogeneous) stress imposed on the array. With regard to the geometry of the array, the 
deformation of the loop is the dominant effect with its radius a deforming as follows: 

a(a) = a (l - x°l°o)- ( 16 ) 

As a result, the self- inductance of the loop L(a) = n aF(a) (with F(a) being a geometry 
dependent factor) will change accordingly: 

L(a)=L (l- Xg cr/a ). (17) 

The relationship between the coefficients x an d Xg is given by 

X 9 = (l + a B g )x (18) 

where 5 9 = ^(f) ao . 

It is also reasonable to assume that in addition to the critical current, the external stress 
will modify the resistance of the contact: 

= ^& - < 19 > 

as well as capacitance (due to the change in the distance between the superconductors): 

C{a) = 1 - C x l/ ao ~ + Xa/ao) - (20) 
To simplify the treatment of the dynamic equations of the array, it is convenient to introduce 
the standard normalization parameters such as the Josephson frequency: 



the analog of the SQUID parameter: 



h ^ ^» (22) 



and the dissipation parameter: 



c 



2nI c (0)C R 2 



(23) 



Combining Eqs. (fUj) and (TTBT) with the stress-induced effects described by Eqs. ( 1T91 and (T2"0"j) 
and using the normalization parameters given by Eqs. fl21l) - (l23l . we can rewrite the equations 
for an array in a rather compact form. Namely, the equations for vertical junctions read: 



1 



+ 



sin ( 



'■j 



+ 



+ 7b- (24) 



Here an overdot denotes the time derivative with respect to the normalized time (inverse 
Josephson frequency), and the bias current is normalized to the critical current without 
stress, 7 6 = I b /I c (0). 

The equations for the horizontal junctions will have the same structure safe for the explicit 
bias related terms: 

1 Vfc e -X°/«o ., 



1 / 

1 - W °o 



+ 



ft 



c 



^ + e W<TO sin^. 



1 



J i-l,3 



,+ 



'i+lj-1 



(25) 



Finally, Eqs. (f2^|) and (1251) should be complemented with the appropriate boundary condi- 
tions 20J which will include the normalized contribution of the external flux through the 
plaquette area rj = 2ir^-. 

It is interesting to notice that Eqs. (1241) and (1251) will coincide with their stress-free coun- 
terparts if we introduce the stress-dependent re-normalization of the parameters: 



(26) 



p L = f3 L (l- Xg a/a )e^, 



(27) 
(28) 



f] = r)(l - 2x(r/a ), 
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(29) 



lb = 7& e 



(30) 



C. Numerical results and discussion 

Turning to the discussion of the obtained numerical simulation results, it should be 
stressed that the main problem in dealing with an array is that the total current through the 
junction should be retrieved by solving self-consistently the array equations in the presence 
of screening currents. Recall [l^] that the Josephson energy of a single junction for an 
arbitrary current / through the contact reads: 



The important consequence of Eq. (l3T]) is that if no current flows in the array's junction, 
such junction will not contribute to the TEC (simply because a junction disconnected from 
the current generator will not contribute to the energy of the system). 
Below we sketch the main steps of the numerical procedure used to simulate the stress- 
induced effects in the array: 

1) a bias point I b is selected for the whole array; 

2) the parameters of the array (screening, Josephson frequency, dissipation, etc) are se- 
lected and modified according to the intensity of the applied stress a; 

3) the array equations are simulated to retrieve the static configuration of the phase 
differences for the parameters selected in step 2; 

4) the total current flowing through the individual junctions is retrieved as: 



5) the energy dependence upon stress is numerically estimated using the value of the 




(31) 



'V,h 



I c sin a ; 



(32) 



total current Ify (which is not necessarily identical for all junctions) found in step 4 
via Eq.flEJ; 
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6) the array energy Ef is obtained by summing up the contributions of all junctions with 



the above-found phase differences 



7) the stress-modified screening currents Iij(T, H, a) are computed by means of Eq.(fl5 
and inserted into the magnetic energy of the array E^ = j^T.ij^I- 



rs \2. 



8) the resulting strain field and TE coefficient of the array are computed using numerical 
derivatives based on the finite differences: 



V 



^[Ei + Ej 



Act 



a(T, H) 



Ae 
AT 



(33) 



(34) 



The numerical simulation results show that the overall behavior of the strain field and TE 
coefficient in the array is qualitatively similar to the behavior of the single contact. In Fig. H] 
we have simulated the behavior of both the small junction and the array as a function of 
the field across the barrier of the individual junctions in the presence of bias and screening 
currents. As is seen, the dependence of a(T, f) is very weak up to / ~ 0.5, showing a strong 
decrease of about 50% when the frustration approaches / = 1. A much more profound 
change is obtained by varying the temperature for the fixed value of applied magnetic field. 
Fig. [5] depicts the temperature behavior of a(T,f) (on semi-log scale) for different field 
configurations which include barrier field / frustrating a single junction and the flux across 
the void of the network r\ frustrating the whole array. First of all, comparing Fig. [5(a) 
and Fig. [2] we notice that, due to substantial modulation of the Josephson critical current 
Ic{T,H) given by Eq.([3l), the barrier field / has similar effects on the TE coefficient of 
both the array and the single contact including temperature oscillations. However, finite 
screening effects in the array result in the appearance of oscillations at higher values of the 
frustration / (in comparison with a single contact). On the other hand, Fig. [5(b-d) represent 



the influence of the external field across the void r\ on the evolution of a(T, / Wreca 



1 that 



in the absence of stress this field produces a well-defined magnetic pattern [18|, [19|, |20| ) . 
As is seen, in comparison with a field- free configuration (shown in Fig. [5(a)), the presence 
of external field rj substantially reduces the magnitude of the TE coefficient of the array. 
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FIG. 4: (Color online) Numerical simulation results for an array 5x5 (black solid line) and a small 
single contact (red dashed line). The dependence of the normalized TEC a(T, f)Tc/\e(0, 0)| on the 
frustration parameter / (applied magnetic field H across the barrier) for the reduced temperature 
T/Tc = 0.95. The parameters used for the simulations: rj = 0, (3 = 0.1, i(0)/£ = 1, £/Al = 0.02, 
Pl = 10, 76 = 0.95, and Xg = X = 0.01. 

Besides, with r\ increasing, the onset of temperature oscillations markedly shifts closer to 
T c . 

IV. CONCLUSION 

We have investigated the influence of a homogeneous mechanical stress on a small single 
Josephson junction and on a plaquette (array of 5 x 5 junctions). We have shown how 
the stress-induced modulation of the parameters describing the junctions (as well as the 
connecting circuits) produces such an interesting phenomenon as a thermal expansion (TE) 
in a single contact and two-dimensional array (plaquette). We also studied the variation 
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FIG. 5: (Color online) Numerical simulation results for an array 5x5. The influence of the flux 
across the void of the network rj frustrating the whole array on the temperature dependence of the 
normalized TEC a(T, /)Tc/|e(0, 0)| for different values of the barrier field / frustrating a single 
junction for 7& = 0.5 and the rest of parameters same as in Fig. 4. 
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of the TE coefficient with an external magnetic field and temperature. In particular, near 
Tc (due to some tremendous increase of the effective "sandwich" thickness of the contact) 
the field-induced TE coefficient of a small junction exhibits clear temperature oscillations 
scaled with the number of flux quanta crossing the contact area. Our numerical simulations 
revealed that these oscillations may actually still survive in an array if the applied field is 
strong enough to compensate for finite screening induced self-field effects. And finally, it 
is important to emphasize that our analysis refers to regular arrays with square geometry 
(similar to already existing artificially prepared arrays However, we can argue that 

the predicted here effects should manifest themselves also in granular superconductors which 
exhibit quite pronounced stress dependent behavior upon mechanical loading 9), l^ . 
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